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ABSTRACT 

A  16  spark-gap  camera  was  used  to  record  the  dynamic  photoelastic  patterns 
of  ten  centrally  cracked,  Homalite-100  specimens  which  fractured  under  ten 
biaxial  stress  ratios  ranging  from  3.7  to  0.  The  dynamic  photoelastic 
patterns  of  curved  cracks  were  used  to  verify  the  previously  developed  dynamic 
crack  curving  criterion.  Cracks,  which  immediately  curved  upon  propagation 
in  three  specimens  under  high  biaxial  loadings,  were  used  to  verify  the  static 
counterpart  of  the  dynamic  crack  curving  criterion.  A  previously  developed 
dynamic  crack  branching  criterion  was  verified  by  the  dynamic  photoelastic 
results  involving  cracks  which  eventually  branched  under  low  biaxial  loadings. 


INTRODUCTION 


Dynamic  crack  curving  and  branching  represent  one  aspect  of  dynamic  crack 
propagation  behaviors  which  has  been  the  subject  of  numerous  theoretical  Cl-51 
and  experimental  C6-9]  investigations.  The  biaxial  load  effects  on  static  and 
quasistatic  fracture  have  also  been  the  subject  of  theoretical  and 
experimental  studies  CIO-121.  The  purpose  of  this  paper  is  to  report  our 
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studies  on  the  combination  of  these  two,  namely,  dynamic  crack  curving  and 
branching  under  various  biaxial  loading  conditions. 

The  crack  curving  criteria  C33  referred  to  in  this  paper  is  a 
micro-mechanic  model  of  continuous  micro-flaw  growth  and  coalescence  in  the 
vicinity  of  the  propagating  crack  tip.  A  similar  model  based  on  qualitative 
evidences  has  also  been  proposed  C133.  The  model  assumes  that  the  crack  kinks 
or  bifurcates,  when  an  off-axis  micro-flaw  connects  with  the  crack  tip.  It  is 
the  dynamic  extension  of  the  crack  curving  criterion  proposed  by  Streit  and 
Finnie  C51.  Thi3  dynamic  crack  curving  criterion  has  been  used  to  predict  the 
crack  kinking  angle  of  a  propagating  crack  under  pure  mode  I  as  well  as  mixed 
mode  conditions  C33.  The  crack  branching  criterion  C43  on  the  other  hand, 
requires  a  critical  stress  intensity  factor  to  trigger  crack  branching  and  a 
crack  curving  criterion  for  predicting  the  crack  branching  angle.  The 
objective  of  this  paper  is  to  provide  further  evidences  in  support  of  the 
dynamic  crack  curving  and  branching  criteria  under  biaxial  loadings. 

DYNAMIC  CRACK  CURVING  CRITERION 

The  authors  have  used  in  the  past,  both  the  maximum  circumferential  stress 
criterion  C 33  and  the  minimum  strain  energy  density  criterion  C143  to  predict 
crack  curving.  The  crack  curving  angles  predicted  by  using  either  of  the  two 
fracture  criteria  are  nearly  identical  for  smaller  values  of  mode  II  stress 
intensity  factor,  Kir  and  thus  only  the  circumferential  stress  criterion  will 
be  used  in  this  paper.  The  near  field  circumferential  stress,  o00,  of  a 
mixed-mode  crack-tip  propagating  at  constant  velocity  is  a  function  of  the 
modes  I  and  II  dynamic  stress  intensity  factors,  K^  and  Kjj,  and  the  remote 
stress  or  the  non-singular  stress  acting  in  the  direction  of  crack 

is  also  a  function  of  the  crack  velocity,  c^. 


propagation. 


dilatational  and  distortional  stress  wave  velocities,  c^  and  c 2,  along  with 

the  other  higher  order  terms  (HOT)  in  the  crack  tip  stress  field.  Crack 

curving  occurs  when  Oqq  reaches  a  maximum  value  at  an  inclination  of  9  and  a 

distance  of  r  away  from  the  propagating  crack.  This  extremum  condition  for 

yields  a  functional  relation  between  the  r  and  9  values.  By  setting  9  =  0  in 

this  relation  under  pure  mode  I  loading,  i.e.  Kjj  =  0,  a  characteristic  radial 

distance,  rQ,  for  self-similar  crack  extension  is  obtained  in  terms  of  Kj, 

o  ,  c,  c.  and  c-,  and  the  higher  order  terms,  HOT. 
ox  1  2 

This  crack  instability  criterion,  which  i3  a  dynamic  extension  to  static 
crack  directional  instability  criterion  of  Streit  and  Finnie  C53,  assumes  that 
the  crack  will  propagate  straight  when  the  above  characteristic  distance  of  rQ 
is  larger  than  a  critical  material  parameter  of  r  .  When  rQ  is  less  than  rc, 
the  crack  suddenly  becomes  unstable  and  will  veer  off  to  the  9c  direction 
which  is  determined  from  the  extremum  condition.  Detailed  studies  on  the 
variations  in  thi3  crack  kinking  angle  with  respect  to  variations  in  fracture 
parameters  are  given  in  References  [3,  143. 


DYNAMIC  CRACK  BRANCHING  CRITERION 

The  dynamic  crack  branching  criterion  requires  as  a  necessary  condition, 
sufficient  released  strain  energy  for  generating  simultaneously,  multiple 
cracks.  Thus,  a  critical  3tress  intensity  factor,  KIb,  is  implicated  in  this 
necessary  condition. 

In  order  to  propagate  the  multiple  crack3  simultaneously,  these  cracks 
must  branch  away  from  the  original  3elf-similar  crack  propagation  path.  Thus, 
crack  curving  is  introduced  as  a  sufficiency  condition.  The  crack  branching 
criterion  can  then  be  summarized  as: 

Kj  >  Kjb  Necessary  condition 


Sufficient  condition 


The  crack  kinking  angle  determined  from  the  above  sufficiency  condition  is  one 
half  of  the  crack  branching  angle. 


FRACTURE  UNDER  BIAXIAL  LOADING 

Biaxial  loading  of  fracture  specimens  have  been  reported  to  influence, 
among  others,  the  apparent  fracture  toughness,  fatigue  crack  propagation  rate, 
plastic  zone  size  and  the  direction  of  crack  propagation.  Of  particular 
interest  to  this  investigation  is  the  predicted  direction  of  the  initial  crack 
extension  of  a  centrally  cracked  plate  under  biaxial  loading  in  Reference 
C123.  Using  the  maximum  circumferential  stress  criterion,  Liebowitz  et  al 
showed  that  the  directional  instability  of  a  central  crack  occurs  when 

r  "aox  >  0.1493  (B-l)  when  B  >  1.0 

0  K 

h 

where  rQ  is  a  small  but  unspecified  characteristic  distance  from  the  crack  tip 

and  B  is  the  ratio  of  the  biaxial  loads,  i.e.,  B  =  F  /F  .  Since  this  crack 

x  y 

kinking  criterion  corresponds  to  the  static  counterpart  of  the  authors' 
dynamic  crack  curving  criterion,  both  criteria  will  be  used  to  evaluate  the 
crack  kinking  data  obtained  in  fracturing  centrally  cracked  plates  with  blunt 
notches  and  under  biaxial  loadings. 


EXPERIMENTAL  AND  DATA  REDUCTION  PROCEDURES 

Dynamic  photoelasticity  was  used  to  determine  the  dynamic  stress  intensity 
factors  in  a  fracturing  Homalite-100  specimen  which  was  loaded  in  a  biaxial 
testing  machine.  The  transient  isochromatics  surrounding  the  propagating 
cracks  were  recorded  by  a  sixteen  spark-gap  Cranz-Schardin  camera.  The  spark 
13  sensed  by  an  EG&G  Lite-mike  and  recorded  on  a  Bioraation  2805M  wave  form 
recorder . 
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The  biaxial  testing  machine  consists  of  two  sets  of  two  oppositely  mounted 
hydraulic  cylinders,  which  allows  the  center  of  the  specimen  to  remain 
stationary  as  shown  in  Figure  1.  Each  set  has  its  own  dented  control  valve, 
accumulator,  handpump  and  pilot  operated  check  valve.  The  two  hydraulic 
systems  can  be  operated  independently,  in  tension  and  in  compression,  as  well 
as  simultaneously  while  maintaining  a  fixed  load  ratio.  The  two  loads,  which 
are  measured  by  two  load  cells  and  through  the  necessary  strain  gage 
conditioners  and  amplifiers,  are  recorded  both  on  a  Brush  recorder  and  on  a 
Biomation  2805M  wave  form  recorder. 

Figure  2  shows  a  16.4  mm  thick,  fractured  Homalite-100  cruciform  specimen 
with  a  central  starter  crack  of  50  ram  length  in  a  test  area  of  approximately 
254  x  254  mm.  The  elastodynamic  properties  of  Homalite-100  were  determined  by 
a  split  Hopkinson  bar  tester  which  yielded  modulus  of  elasticity,  Poisson 
ratio  and  stress  optical  constant  of  4.27  GPa,  0.357  and  32.8  MPa-mm/fr, 
respectively  at  a  strain  rate  of  about  2  x  10  1/sec. 

The  least  square  procedure  C15,161  was  used  to  determine  the  fracture 
parameters,  i.e.  the  stress  intensity  factors,  K-j-  and  Kjj,  where  applicable, 
the  remote  stress  component,  jqx,  and  two  additional  higher  order  terms  (HOT) 
of  Mode  I  and  Model  II,  A^  and  B^,  respectively  from  the  transient 
isochromatics.  Convergence  of  thi  procedure  was  verified  for  each 
computation  by  the  goodness  of  fit  between  the  experimentally  obtained 
isochromatics  and  those  reconstructed  from  computed  dynamic  fracture 
parameters . 

RESULTS 

A  total  of  ten  specimens  with  different  biaxial  load  ratios,  as  shown  in 
Figure  3,  were  tested.  The  ratios  of  applied  load,  which  are  parallel  and 
normal  to  the  central  crack,  ranged  from  B  -  3.7  -  0. 
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The  45  degree  crack  in  Specimen  No.  DH073084,  which  was  loaded  to  failure 
under  nearly  equal  biaxial  load,  i.e.  B  =  1.07,  propagated  in  a  self  similar 
straight  line,  thus  indicating  that  the  biaxial  testing  machine  generated  a 
uniform  stress  field.  On  the  other  hand,  the  45  degree  crack  in  Specimen  No. 
DH080284,  which  failed  at  a  biaxial  load  ratio  of  B  =  1.53,  immediately  curved 
upon  propagation  as  expected.  Dynamic  crack  curving,  which  occurred  in 
specimen  under  low  biaxial  load  ratios  of  B  =  1.6,  1.33  and  1.16  in  Specimen 
Nos.  DH072084,  DH062284  and  DH062584,  respectively,  were  also  used  to  evaluate 
the  dynamic  crack  curving  criterion.  Cracks  emanating  from  the  horizontal 
cracks  in  Specimens  Nos.  DH071684,  DH071784  and  DH071984  under  biaxial  load 
ratios  of  B  =  3.7,  3.06  and  2.6,  respectively  immediately  curved  upon 

propagation  and  before  the  camera  could  be  triggered.  These  results  were  used 
to  evaluate  the  authors  C31  and  Liebowitz  et  al  C12J  static  crack  curving 
criteria.  Figure  4  shows  Specimen  No.  DH071784  in  which  the  crack  branched 
immediately  upon  propagation  rrom  two  sharp  starter  crack  tips.  The  crack 
branching  data  which  occurred  after  crack  curving  in  three  specimens,  which 
includes  the  above  Specimen  No.  DH071784,  were  used  to  evaluate  the  crack 
branching  criterion.  Selected  results  are  discussed  in  the  following. 

Figure  5a  shows  the  dynami;  photoelastic  fringe  patterns  associated  with 

Frame  No.  ?  of  Specimen  No.  BH0/2064.  Figure  5b  shows  the  data  points  used  in 

evaluatinq  KT,  KTT,  a  ,  A,  and  E-,,  as  well  as  the  reconstructed  isochromatics 
I  II  ox  3  3 

usina  these  crack  tip  parameters.  Figures  6a  and  6b  shows  the  variations  in 

Kt ,  KTr  and  1  in  the  right  and  left  cracks  of  Specimen  No.  DH072084. 

I  II  ox  ^  r 

Figure  7  shows  the  crack  path  which  resulted  in  a  crack  branching  in 
Specimen  No.  DH072084.  Also  shown  in  Figure  7  is  the  corresponding  measured 
and  predicted  crack  curving  angles  along  the  crack  path.  The  latter  was 
computed  by  substituting  the  resultant  K^.,  and  and  the  two  higher 
order  terms,  A  and  E,,  into  the  crack  curving  criterion. 

j  J 


Figure  8  shows  two  frames  of  the  dynamic  isochromatics  prior  to  and  after 
crack  branching  in  Specimen  MD012884.  Figures  9a  and  9b  show  the 
variations  in  Kj,  and  oqx  before  crack  branching  of  the  left  and  right  cracks. 
Extrapolated  and  oqx  values  were  then  used  to  determine  the  dynamic  crack 
branching  stress  intensity  factor  as  well  as  the  branching  angle. 

Table  1  shows  the  measured  and  computed  crack  curving  angles  immediately 
a. ter  crack  propagation  in  the  five  specimens  3hown  in  Figure  3.  Both  the 
authors'  [31  and  Liebowitz  et  al  [121  procedures  were  used  to  estimate  the 
initial  crack  curving  angles.  Good  agreements  between  al]  theoretical  and 
experimental  results  are  noted. 

Table  2  3hows  the  measured  crack  branching  stre33  intensity  factor  as  well 
as  the  measured  and  estimated  crack  branching  angles  in  four  specimens.  At 
the  lower  biaxial  ratios,  B  <  1.0,  the  experimentally  observed  crack  branching 
angles  are  within  24  -  44  degrees  as  reported  by  others  [4,8,16-181  under 
uniaxial  loading.  The  crack  branching  angles,  however,  increased  with 
increase  in  the  biaxiality  ratios.  It  is  also  interesting  to  note  that  a 
meager  13%  increase  in  the  branching  stress  intensity  factor,  resulted  in 

multiple  crack  branching. 

The  two  Tables,  I  and  II,  show  that  the  curving  and  branching  angles  can 
be  predicted  with  reasonable  accuracy.  Regardless  of  the  biaxial  loading  ratio 
and  the  branching  crack  length,  the  branching  stress  intensity  factor  of  2.0 
MPa1 m  is  obtained  for  Homalite-100  and  is  in  good  agreement  with  the  results 
reported  in  Reference  [43. 

CONCLUSIONS 

1.  Dynamic  crack  curving  branching  criterion  proposed  by  the  authors 

successfully  predicted  the  crack  curving  and  crack  branching  angle  under 
;  extremes  observed  in  biaxial  loading. 


8 


2.  A  crack  branching  stress  intensity  factor  of  K 


2.0  MPa^  was 


lb 

determined  for  this  Homalite-100  sheet  and  is  independent  of  biaxial 
loading  ratio. 
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TABLE  I 

CRACK  CURVING  ANGLE  (in  degrees) 


TEST 

SPECIMEN 

NO. 

Biaxial 

Load 

Ratio  B 

Measured 

0 

m 

Predicted 

9t 

Ramulu  et  al 

Predicted 

et 

Liebowitz  et  al 

DH071684 

3.7 

76 

71.7 

72.5 

DH071784 

3.05 

67 

66.2 

65.7 

DH071984 

t* .  6 

54 

52 

56.5 

DH072084 

1.6 

5 

0 

0 

DH080284 

1.53 

30 

32.5 

29 

LI 


A 
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TABLE  2 

SUMMARY  OF  CRACK  BRANCHING  DATA  UNDER  BIAXIAL  LOADING 


CRACK  BRANCHING  ANGLE  IN  DEGREES 


SPECIMEN 

NO. 

Biaxial 

Load 

Ratio  B 

( MPa  *  in ) 

KIb 

Measured 

r» 

m 

Estimated 

e 

MD012284R 

0.0 

1.87 

28 

40 

MD012284L 

1.98 

24 

22 

MD0L28S4R 

0.33 

1.98 

44 

44 

MD012884L 

2.05 

40.5 

40 

DH072084R 

1.6 

2.0 

49 

- 

DH072084L** 

2.1 

DH071984R 

2.6 

1.98 

50 

- 

DH071984L** 

2.12 

■'H071784R* 

3.05 

- 

- 

DH071 7  ,;4L** 

2.12 

i.verace 

2.02 

L:  Left  Crack 
R:  Right  Crack 
*  No  Branchina 
**  Multiple  Branching 
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